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Abstract 

We investigate a particular phase transition between two different tunneling regimes, direct and 
injection (Fowler- Nordheim) , experimentally observed in the current- volt age characteristics of the 
light receptor bacteriorhodopsin (bR). Here, the sharp increase of the current above about 3 V 
is theoretically interpreted as the cross-over between the direct and injection sequential-tunneling 
regimes. Theory also predicts a very special behaviour for the associated current fluctuations 
around steady state. We find the remarkable result that in a large range of bias around the 
transition between the two tunneling regimes, the probability density functions can be traced back 
to the generalization of the Gumbel distribution. This non-Gaussian distribution is the universal 
standard to describe fluctuations under extreme conditions. 

PACS numbers: 05.70.Fh,05.06.Gg 
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INTRODUCTION 



Non Gaussian distributions (NGDs) evidence deviations from the central limit theorem, 
ubiquitous under thermal equilibrium conditions. Among the NGDs, increasing attention 
has been recently devoted to the case of the generalized Gumbel distribution G(a), related 
to the statistics of extreme events. Since the pioneer papers p~H3] , Gumbel-like distributions 
have been evidenced in a wide series of extreme events, ranging from physical or sociological 
and environment scenarios [IH9]. The possibility that the generalized Gumbel distribution 
can be associated with the identification of some universal behaviour of a wide class of 
extreme events is emerging as an intriguing issue [31 EE CO E]- 

The aim of this paper is to present evidence of the G(l) distribution in the current 
fluctuations of a monolayer of bacteriorhodopsin (bR), a protein present in the Archean 
organism Halobacterium salinarum[10\ and acting as a light receptor. 

In the framework of a microscopic interpretation of electrical transport through this 
protein, a model, called INPA (impedance network protein analogue) [TT| [T2] . reproduces 
the current- voltage (I-V) characteristics in agreement with experiments [13]; on the other 
hand, it also produces the current fluctuations on which no data are presently available. 
A set of I-V characteristics were carried out on bR monolayers, making use of an atomic 
force microscope (AFM) technique. Measurements were performed at the nano-scale length 
and evidenced a cross-over between a direct (DT) and an injection or Fowler-Nordheim 
(FN) regime of the associated tunneling current [13]. The aforementioned INPA succeeded 
in providing a microscopic interpretation of experiments, by using an interpolation scheme 
that accounts for a continuous transition between the two tunneling regimes. The mechanism 
of charge transfer among the amino-acids constituting the tertiary structure of the protein, 
is taken to be a sequential tunneling. 

Figure [T] shows a sketch of the experimental apparatus. The I-V characteristics, as ob- 
tained by the experiments (see continuous curve) and the theoretical model (see the full 
circles) when the AFM tip just touches the protein films at about 4.6 nm from the bottom 
metal-contact, is also reported. The transport approach uses the Monte Carlo technique and 
allows for the simultaneous detection of the current and its fluctuations around the steady 
state. 

Interestingly enough, we found that the cross-over between DT and FN regimes is much 
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strongly marked by the variance of current fluctuations than by the current variation itself. 
Moreover, the associated NGDs are found to follow a generalized Gumbel distribution [3J [5] ■ 
The content of the paper is the investigation of the properties of the current fluctuations 
and is organized as follows. In the next section we will briefly recall the microscopic model 
used to interpret the experiments [13]; this model predicts NGDs for the current fluctuations 
around steady state. Section 3 will briefly summarize the properties of the scaled Gumbel 
distribution, G(l). Section 4 will report the results of the theoretical predictions and discuss 
the implications of the NGDs of current fluctuations exhibited by the simulations. Major 
conclusions will be reported in Sec. 5. 



THEORY 



The layout of the microscopic INPA approach is briefly summarized in the following. The 
single protein is mapped into an impedance network, by using the C a atom of each amino 
acid as a node and by introducing an elemental impedance (in general an RC circuit) as a 
link between neighbouring nodes. The number of links starting from each node is determined 
by a cut-off interaction radius, say Rc, and by the network topology (not regular). In the 
present analysis, we choose Rc = 6 A, a value that optimizes the native to activated state 
resolution of the protein [TT]. The value of the elemental impedance (a simple resistance in 
this case) depends on the distance between amino-acids as: r^j = pkj/Aij, where p is the 
resistivity. The value of resistivity is here taken to be the same for all the links, and depends 
on the voltage drop between the corresponding nodes, as detailed below. The pedices i,j 
refer to the amino-acids between which the link is stretched, 1^ is the distance between the 
labeled amino-acids, taken as point-like centers, and Aij is the cross-sectional area shared 
by the labeled amino-acids: Aij = vr ( R c — ^j/4j [Hj. 

To take into account the superlinear features of the I-V characteristics [13], the link 
resistivity, p, is chosen to depend on the voltage drop between nodes as: 

Pmax eV < $ 

P(V) = (1) 

where Pmax is the resistivity value used to fit the I-V characteristic at the lowest voltages, 
Pmin "C Pmax is an extremely low series resistance limiting the current at the highest 



voltages, and $ is the height of the tunneling barrier between nodes. Since charge transfer 
is here interpreted in terms of a sequential tunneling between neighbouring amino-acids, the 
above interpolation formula reflects the different voltage dependence in the prefactor of the 
current expression [15J: / ~ V in the DT regime, and / ~ V 2 in the FN tunneling regime. 

For the transmission probability of the tunneling mechanism we take the expression given 
by Ref. [16J. The microscopic model of electron transport used for describing the experimen- 
tal data, is based on the local possibility of choosing between different tunneling probabilities 
given by: 



V?, = exp 



<r>FN 

Pij = ex P 



2h d L I 



2m($ - -eVil 
'2l ij V2m\ $ /$ 



eV i}j < $ (2) 



eV id > $ ; (3) 



\ h J eVij V 2 

where Vy is the potential drop between the couple of i, j amino-acids and m is the electron 
effective mass, here taken the same of the bare value. 

Figure [2] reports the evolution of the simulated current responses for two different voltage 
values, 1 V corresponding to DT regime, and 5 V corresponding to FN regime. We can 
observe the increasing number of spikes when passing from the low to the high voltage. At 
5 V the number of spikes is so high that it becomes difficult to classify them as "extreme 
events". On the other hand, by lowering the bias value, the number of spikes lowers but 
does not definitely disappear, as shown in the same figure for the case of 1 V. 

Figure ^ reports the variance of current fluctuations corresponding to the I-V char- 
acteristic of Fig. Q. Here we notice a rather abrupt increase of the variance of current 
fluctuations in concomitance with the cross-over region of the I-V characteristic. The giant 
increase, for about five orders in the magnitude of current variance, is associated with the 
opening of low resistance paths between contacts, originated by the establishing of the FN 
tunneling regime. 

Figure [I] reports the probability distribution functions (PDFs) of current fluctuations at 
different voltages covering the cross-over region of the two tunneling regimes. The variables 
have been rescaled as suggested in [21 El [TF] to evidence the "universal" behaviour of all the 
curves. In the inset of Fig. Q the PDFs are shown before rescaling. The divergence from 
Gaussian behaviour is not completely surprising, being partially due to the small size of the 
system under consideration [17] . but also to the critical condition of a system in which, for 
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all the considered voltages, two different tunneling regimes are in competition. On the other 
hand, this kind of systems, and their fluctuations in particular, have received in the last 
years a relevant attention, mainly due to the possibility to trace back them to an unique 
root. The amount of results of recent past [2H9j suggests this behaviour should be described 
by a generalized Gumbel distribution, one of the reference distributions for the statistical 
modeling of extreme events [18] . 

THE SCALED GUMBEL DISTRIBUTION 

The results reported in the inset of Fig. Q are found to exhibit a skew shape, resem- 
bling the popular, almost ubiquitous PDF studied by Bramwell-Holdsworth-Pintor (BHP) 
[2]. This distribution was found in very different contests going from fluidynamics, to self- 
organized critical systems, and to resistance fluctuations Therefore, we analyze 
current fluctuations to understand whether BHP can be, also in the present case, an ade- 
quate PDF. To this purpose, data are scaled by using the mean value and the variance in 
the conventional way [2]: 

z -> ^— <z> ^> = Xj n(x) ->■ an = G{a, x) (4) 

a 

with < x >, cr assuming the standard meaning of the PDF, the mean value of 

the variable and the standard deviation, respectively, and with z = ln(I). The numerical 
parameter a is related to the specific system conditions [5]. 

The choice of the logarithm of current instead of the current itself is suggested by the 
slow dependence of the former on the applied voltage [2j [19]. As illustrated by Fig. ([5]), 
after rescaling, PDFs related to different bias values (1 4- 5 V) collapse in a single one, say 
G(l,x) which, on the other hand, is not the BHP but the generalized Gumbel function of 
order 1: 

G(l,x) = 6(a)e-^ x+ ^- e - Wa)x+ ' l) (5) 

with a = 1, 8 2 (a) the trigamma function and 7 the Euler constant. This PDF has been also 
called the scaled Gumbel distribution function, and has been found of interest in describing 
the fluctuations of the conditional galaxy density [19] and is a special case of the PDF 
initially introduced by Chisel- Fortin-Holdsworth [7] (CFH). For the sake of completeness, 
the CFH distribution recovers the BHP when a = tt/2. The CFH distribution has been 
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found in many systems in which at least one phase transition appears [51 [8] and usually, as 
made for BHP, its origin was related to the existence of a correlation length larger than the 
system size. On the other hand, it is the natural generalization of the Gumbel distribution 
of extreme events pQ; therefore the question of how the two interpretations can converge 
remains an open one. The question is discussed in Ref. [5], and the CHF distribution with 
integer values of a is described as the distribution of the a-th larger /smaller value of the 
variable, here identified in the high and low resistivity values. The analogy with natural 
(for example the climate) extreme events is evident: they are considered extreme or exotic 
under standard environmental conditions, otherwise when environmental conditions change 
in a significative way, they becomes the most probable response [20J. 

RESULTS AND DISCUSSION 

The INPA model was used in the " large tip" version that involves several amino acids in 
the contact between the protein and the tip. In such a way, the real extent of the AFM tip 
(transversal size of about 200 nm) can be taken into account. The single protein results fit 
the sample response (see Fig. 1) when the resistivities to be associated with the links are 
taken as Pmax = 8 x 10 13 Vt A, p min = 4 x 10 6 £7 A, which correspond to a rough estimation 
of A" w 10 9 trimers in an area of about 2 x 10 n nm 2 [12]. The energy-barrier height is 
$ = 219 meV. We remark that after the proper scaling, the PDFs of current fluctuations 
shown in Figs. (2) and (3) converge to the single G(l,x) curves, as reported in Fig. Q 

The sharp cross-over between the DT and FN tunneling regimes admits for a phase 
transition interpretation, with DT and FN being the two phases. To this purpose, Fig. (J5]) 
reports the tunneling transparency used to model the sequential tunneling mechanism (see 
Eqs. (2) and (3) ), here called probability distribution, vs the applied voltage. The cross- 
over between the DT and FN regimes occurs in the region where the two curves (dashed and 
continuous lines) cross each other, respectively at the two voltages: the critical voltage, V c , 
and, with the same meaning of the Ginzburg temperature, [21], the Ginzburg voltage, Vc- 
For voltages in the region V < V c , the resistivity of each link can take only the two values 
Pmax and p m in, respectively. This is the case of a squared shape for the tunneling barrier. 
For voltages V > V c , which corresponds to a triangular shape of the tunneling barrier, the 
single value of Pmax smooths to a continuous range of p(V). Accordingly, we conjecture 
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that the region of bias inside the values V c Vq corresponds to the consolidation of the FN 
tunneling regime, the new phase for the charge transport through the physical system. In 
other words, by approaching the voltage region V c -j- Vq, for V it j in the DT regime, each 
link will choose a resistivity values between pmax and p m in- By contrast, for Vij in the FN 
regime, each link will choose a resistivity values between p(V) and p m in- For voltages above 
the value of Vq, the FN phase will take place for all the links. 

CONCLUSION 

The paper reports a microscopic model of current- volt age and current fluctuations in a two 
terminal sample of nanometric width where the active region is provided by a monolayer 
of bacteriorhodopsin. In particular, charge transport has been analyzed in a very wide 
range of voltages up to 7 V. Here current response exhibits a strong superlinear behavior at 
increasing voltage, typical of a transport controlled by a tunneling mechanism. A unifying 
microscopic description, the so-called INPA model, is able to interpret the experiments on 
the basis of two tunneling regimes, a direct one (DT) at low bias and an injection one (FN) 
at high bias. Both tunneling regimes are assumed of sequential type and they occur between 
neighbouring amino-acids. Within this model, we obtain good agreement between theory 
and I-V experimental characteristics. More interesting, current fluctuations are found to 
exhibit a sharp increase for about five order of magnitude in correspondence of the cross- 
over between the DT and FN tunneling regimes. The associated probability distribution 
function is found to follow a universal behaviour that is characterized by a generalized 
Gumbel distribution, which is typical of systems close to a phase transition. In the present 
case, the phase transition is associated with the cross-over between two tunneling regimes 
each of them taken as a single phase, and a kind of phase diagram is individuated in the 
voltage dependence of the tunneling transparency that has been used to interpolate the two 
tunneling regimes. 
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FIG. 1: A schematic illustration of the experimental apparatus used for the I-V characterization. 
The experiments (continuous curve) and the theoretical results (full circles) are reported in the 
center of the figure. 
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FIG. 2: Evolution (in units of iterations) of the instantaneous current responses for a low bias of 
1 V, when direct tunneling regime occurs, and for a high bias of 5 V, when injection tunneling 
regime occurs. 
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FIG. 3: Variance of current fluctuations normalized to the value of the steady state current vs 
applied voltage. The continuous curve is a guide to the eyes. 
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FIG. 4: Scaled PDFs of current fluctuations for different bias values. Symbols refer to simulations, 
continuous line to the scaled Gumbel distribution of order 1 (see Eq. 5 in text). The inset reports 
the PDFs at the given bias before applying the scaling. 
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FIG. 5: Tunneling transparency given by Eqs. (2) and (3), here interpreted as probability distri- 
bution of the the resistivities pertaining to the two tunneling phases, vs normalized voltage. 
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